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Abstract 

A small perturbation analysis, in the long wavelength regime, is used to 
obtain the downstream boundary condition for the pressure for the flow over a 
flat plate. The methodology is extendable to other geometries. Numerical 
results for high Reynolds number laminar flows show great improvement in 
convergence rate to steady state as well as in the quality of the results. 
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1 . INTRODUCTION 


Since the early 1970's there has been growing attention paid to the 
numerical solution of the full time- dependent, viscous, compressible Navier— 
Stokes equations of fluid mechanics; see inter-alia [1], [4]. A vexing and 
still open issue is the formulation of "downstream" or "outflow" boundary 
conditions. One common strategy for flows evolving to steady state has been 
to assume that the particular geometry considered allows the imposition of the 
infinity-downstream steady state conditions at the relevant boundary of the 
computational domain. For example, for flow in a pipe, the downstream 
condition might be taken from the fully-developed, steady-state parabolic 
profile. For a semi-infinite flat plate at zero angle of attack, it is 
usually assumed that "sufficiently far" downstream, the free stream conditions 
predicted by boundary layer theory are sufficient for stating the 
corresponding boundary conditions. 

Numerical experiments (see e.g., [5-6]) with time-consistent codes have 
shown that taking P = P OT for the downstream boundary conditions leads to 
very slow convergence to steady state. This phenomenon is aggravated with 
increased Reynolds numbers, because the spatial resolution of the boundary 
layer then leads, via the von Neumann stability analysis, to very small time 
steps. It can be seen from the results of Rudy and Strikwerda [5] that the 
retardation of the convergence to steady state is caused mainly by the slow 
decay of relatively long wavelength disturbances. They tried to account for 
this by proposing a boundary condition on the pressure of the form 
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where p is the density, c = (yp/p) ^ is the speed of sound, and a is a 
free parameter which is "tuned" so as to optimize the convergence rate to 
steady state. Note, that in any case at steady state, the downstream boundary 
pressure has the free-stream value. Their treatment greatly accelerated the 
convergence rate; its derivation, however, is basically "hyperbolic." 

In this paper we propose to derive downstream boundary conditions using a 
different methodology. The basic idea is to perturb the Navier-Stokes 
equations around some approximation (usually a fairly crude one) to the steady 
state solution in the downstream region. The linearized partial differential 
equations for the perturbed quantities are then attacked using a modal form 
for their solution. The resulting ordinary differential equations present an 
eigenvalue problem which is solved assuming very long wavelengths. This 
eigenvalue problem yields the decay rate for long waves and their phase 
velocity. With these two quantities, as functions of the Mach number and the 
Reynolds number, one can use the modal form of the perturbed solution to 
obtain the desired boundary condition. 

In section 2 we give the derivation of the eigenvalue problem and the 
boundary condition for the particular case of the flow past a two-dimensional 
flat plate. 

Numerical results are presented in section 3 for the case of subsonic 
Mach number, M ra = 0.4. The corresponding Reynolds numbers per foot 
is 3 x 10 5 . The numerical experiments are carried out for semi-infinite 
plates with the outflow boundary taken at various locations downstream of the 
leading edge. 



2. THEORETICAL DEVELOPMENT 


The geometry under consideration is shown in the sketch below. 



The two-dimensional compressible Navier-Stokes equations governing the flow 
may be written in conservation form as 
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The various parameters and dependent variables P> u, v, y, X, E, k, T, and 
P are > respectively, the density, the x and y components of velocity, the 
shear and second coefficients of viscosity, the total energy per unit volume, 
the coefficient of heat conduction, the temperature, and pressure. An 
equation of state p = p(e,p) relates the pressure to the density and the 
specific internal energy e, 


e = (E/p) - j ( u 2 + v 2 ). (2.3) 

For example, for ideal gases the equation of state is 

p = (y-1 )p[e + — (u 2 + v 2 )] (2.4) 

where Y is the ratio of specific heats at constant pressure and volume, 
respectively. 

The numerical algorithm used in the next section is based on solving a 
set of difference equations which approximate the full Navier-Stokes equations 
(2.1) — (2.4). However, in considering the perturbation model we make several 
rather crude simplifications. 
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We assume, first of all, that near the outflow boundary the steady state 
is given by a constant pressure p Q , constant density p Q , a vanishing 
vertical velocity Vq = 0, and that the velocity in the x-direction is given 
by 


u 


0 



(6 < y < °°) 
(0 < y < 6) 


(2.5) 


where 6 is the nominal edge of the boundary layer and a = x/p is the 
velocity gradient at the wall (y = 0). The shearing stress, x, is, of course, 
a function of x but will not be treated as such in the development. 

A second simplification is to assume (for the mean steady flow) that the 
total enthalpy, Hq, is constant, i.e.: 
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Equation (2.6) is equivalent to the assumption that we are perturbing around a 
steady flow of Prandtl number unity. 

The perturbed field quantities are defined as departures from the assumed 
steady state described above: 


u = u - u. 


P = P “ V f 


P = P “ P, 


(2.7) 
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Thus, p' is the pressure perturbation, etc. Next, we linearize the Navier- 
Stokes equations in the boundary layer region (0< y< 6) by substituting 
(2.7) into (2.1) and retaining only first-order quantities in the primed 
variables. The resulting set of equations may be written as follows: 


!t“ + ay fir + p 0^fx“ + It*) = °* (continuity) (2.8) 
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(x-momentum) (2.9) 


(y-momentum) , (2.10) 


and instead of the linearized energy equation we have the linearized constant 


enthalpy relationship: 
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In the above equation is the kinematic viscosity at steady state 

conditions. We have also used the thin layer approximation in (2.9) and 
(2.10). We emphasize again that the perturbation equations (2.8) - (2.11) are 
valid only for 0 < y < 6. Outside this viscous region we shall derive 
another set of equations which then will have to be matched to the above 
equations. First, however, we set the modal analysis. Because (2.11) 
provides an algebraic relationship between p' and p" and u' , it is 
necessary to consider only a three-vector, (p', u" , v") for which we set the 


modal ansatz: 
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e i( 4 »t+bx) 



( 2 . 12 ) 


Next, we normalize and nondiraensionalize the various quantities as follows: 


n = ay/U w , 
w = ty/a , 


G 1 ' F l /p -' G 2 ‘ r 2 /U -- G 3 - V°„. 

6 - Ub/a, t - - 1/R { , 


where is the Reynolds number based on the boundary layer thickness. 

In terms of the dimensionless variables and parameters, after 
substituting ( 2 . 12 ), equations ( 2 . 8 ) - ( 2 . 10 ) become, for 0 < n < 1 : 

i(io + Bn)G + igG 2 + 6 ' = 0 (2.13) 

l(w + 7 h)G 2 + -^2 Gj + G 3 - eG" (2.14) 

l(w + Bn)G 3 +-ijGj - III nG ' - G 2 = eG^. (2.15) 

^ ^oo 

In equation (2,13) - (2.15) the primes designate differentiation with respect 
to the coordinate n; e.g., = dG^/dn, G ^ 1 = d 2 G^/dq 2 . 

Next, we consider the inviscid region, q > 1 where Uq = and the 

flow is described well by the Euler equations (e = 0). Repeating the above 
procedure the acoustic perturbation equations corresponding to (2.13) - (2.15) 


are derived: 
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i(u + B)G 1 + i|3G 2 + G' = 0 (2.16) 

i(u) + £)G 2 + Gj = 0 (2.17) 

' Y M „ 

i((0 + 6)G 3 + G 1 " ^ g; = 0. (2.18) 

ym y 

00 

The solvability conditions for the existence of solutions of the form e rr| 
for the above homogeneous set of differential equations is found, after some 
manipulation, to be 

r 2 = 6 2 - yM 2 (w + £■)(« + 6). (2.19) 

The set of equations (2.13) - (2.15) constitutes a fifth-order system and 
requires five boundary conditions. Two of them are supplied by the no-slip 
condition at the plate: 

G 2 (0) = 0 (2.20) 


G 3 (0) = 0. 


( 2 . 21 ) 


The other three are supplied by requiring the correct "impedance-matching" at 
P = 1 between the system (2.13) - (2.15) and the set (2.16) - (2.18), namely: 


Gj(l) = rG x ( 1 ) 

(2.22) 

G'(l) = rG 2 (l) 

(2.23) 

G'(l) = rG 3 (l), 

(2.24) 
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where r = r(M ; to ,3) is given by (2.19). 

CO 

The system (2.13) - (2.15) with the boundary conditions (2.20) - (2.24) 
constitutes the eigenvalue problem from which, in principle, we could extract 
the dispersion relation oi = w(8) for given Mach and Reynolds numbers 
(M^ and e** 1 ). 

As mentioned in the introductory section, we are interested in very long 
waves, i.e., the small 8. We, therefore, propose expanding the eigenfunc- 
tions Gjl (i = 1,2,3) in the form 

G ± = G< 0) + & g [ 1) + 8 2 G< 2) + (i = 1,2,3) (2.25) 

where 

G i°^ = [ G i(n;“»3,e)] B= o. (2.26) 

We also expand the dimensionless "frequency” w (which for small 0 will 
turn out to be the decay rate) in a similar manner: 

a) = o)q + 0u)^ + ••• (2.27) 

where 

(2.28) 

Using the expansions (2.25) and (2.27) in (2.13) - (2.15), we get the 

following zeroth and first-order problems in 0: 
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(2.29) 

!< 0) - eG< 0) " = 0 

(2.30) 

17 (nc< 0> y - £ c<°>" - 0 

(2.31) 


and 


im Q + G^' = 


- i» 1 g<°> - iG<°> - i nG < o > 


(2.32) 


1U 0 G 2 n + G 3 1} " £G 2 1) " = " + 7j G 2 0) " *^2“ G l 0) (2 ‘ 33) 


*» 0 <4° + ~±2~ G< lj - 2^- - eG^ 1 ^" = - Ku^ + n)G 3 (0) . (2.34) 


Similar sets of equations can be derived for higher powers in g. They 
all have the same left-hand side in the relevant variables, with right-hand 
sides serving as source terms made up of all previous (presumably known) 
approximations. It is seen that u)q, the zeroth-order approximation to the 
decay rate is determined from the "g = 0 problem," equations (2.29) - (2.31), 
without reference to oj^ = (dw/dg)g_g, the phase velocity of the perturba- 
tions. We, therefore, first concentrate our attention on this "g=0 
problem." By some simple manipulations one can extract from the system (2.29) 
- (2.31) a single fourth-order differential equation for G^. For the sake 
of neater notation set G^ - g ^ , (j = 0 , 1 ,...). With this change of 
notation, the fourth-order differential equation for g Q = G^ is found to 
be 
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The boundary conditions for (2.35) are, to zeroth order in B: 


g 0 (0) = g^O) = 0 


(2.36) 


g^D = ±o 


V2 


y-l 
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O3 0 g 0 (D 


(2.37) 


g^'U) = io V 2 ^ V 2a) 0 gJJ(l). (2.38) 

9 

In equations (2.35) - (2.38), a = (y-DM^. For high Reynolds number flow, 

say Rg > 30, e is a small number, and we can ask ourselves how does the 
eigenvalue behave as a function of e as e diminishes and approaches 

zero • 

There are three possibilities: 


0 
00 

constant 

Another way to pose this is to say that u)q ~ e a as e 0 and we ask whether 
a is smaller, greater than, or equal to zero. If a < 0, then we have in 
(2.35) a singular perturbation problem and one can show that, considering the 
boundary conditions (2.36) — (2.38), the only way to match the inner and outer 
expansions is to have the trivial solution. If a > 0, we have a regular 
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perturbation problem, and again the only possible solution is the trivial one 
of gg(h) 0 as e 0. We thus remain with the third possibility of 
u 0 ~ e ^ or sma H £• Let this proportionality constant be s/ujq = Xq. 
Then (2.35) is rewritten as: 



(2.39) 


It is assumed, and will be verified, that for reasonable Mach numbers, say 
below hypersonic speeds, X Q = 0(1). 

We can now state the "3 = 0 problem" more concisely: given o, e, and 

Y (i.e., the Mach number, Reynolds number, and the gas constant) solve the 
differential equation (2.39), subject to boundary conditions (2.36) - (2.38), 
for the eigenvalues Xq and the eigenfunction gQ(n) = This 

eigenvalue problem was solved numerically for a range of values 0 < e < .05, 

0 < o < 20 (for laminar flow this corresponds roughly to Reynolds 

number, 50 < < «> and Mach number, 0 < M^ < 7 for y = 1.4). For this 

range of the parameters it is found that X Q is purely imaginary negative — 
this means that (d q is purely imaginary positive. Thus, to Q represents the 
decay rate of very long waves — to be corrected for finite but long waves by 
the term 8a)j. Figure 1 shows a plot of |1 /Xq| versus a; for 
0 < e < .0053 (10 < ^bo < 00 ), Xq is completely insensitive to changes in the 
Reynolds number and indeed may be taken to depend on the Mach number only. 
As a + 0, I l/Xg | takes on the value tt 2 /4 which can be derived 
analytically. 

It is seen from Figure 1 that, for this range of parameters, |X Q | is 
indeed of order unity, .4 < |X Q | < 5. For R^ as low as 40 (not shown on 
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the graph) the j 1 /X q | curve deviates from the one shown by less than 5%. 
Thus, we conclude that for moderate and high Reynolds numbers the very long 
wave (3 ~ 0) decay rate is basically a universal function of a = (y~1)m£. 
The correction, 3Wj, to the decay is found by considering the first-order 
problem, (2.32) - (2.34), with its attendant boundary conditions. Again, one 
can extract from this system a single fourth-order differential equation for 
gj = After some manipulations the g^-problem may be stated as follows: 
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The boundary conditions, to first-order in 3, are found to be 


g x (0) = g'i'(O) = 0 


g'(l) - ia 
g'"(l) - ia 
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where in (2.40) - (2.46) for given a and s we already know, albeit 
numerically, gg and uig. In deriving the boundary conditions (2.45), (2.46) 
one has to expand the expression (2.19) for r in powers of 0. One then 
finds that the expansion is valid only for 0 /<j>q small compared to unity. 
(Exactly how small cannot be answered at this point without finding the radius 
of convergence for our expansion.) In physical terms, for laminar flow where 
6 ~ /x, this constraint means that our analysis is restricted to modes whose 
wavelength is longer than the computational domain on the plate. 

A solvability condition for the gj-problem is found by multiplying (2.40) 
by the left null vector of the left-hand side operator of (2.40). 
Symbolically, we may write 


<f,R 2 > 
“l = " <f 


(2.47) 


where f is the left null vector to the system, including the boundary 
conditions. Numerical evaluation of (2.47), for the same range of the 
Reynolds and Mach numbers mentioned before, leads to the result that oi^ is 
almost a universal constant in the present problem, namely 

ojj = - .295 ± .002. (2.48) 

Now that we have the decay rate oJq and the phase velocity for long 

waves, we can derive outflow boundary conditions. 

Consider, e.g., the pressure perturbation: 
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p' ~ e iin + bx) G(n) (2.49) 

where G is a linear combination of Gj and G£. Since the numerical 
solution of the finite difference approximation to the Navier-Stokes equations 
is often done in terms of dimensionless temporal and spatial variables, we 
define dimensionless time and x coordinates by 

U t 

e - “Yj - ; K = £ (2.50) 


where the reference length L is the distance from the leading edge to the 
outflow. Using \J> = au, b = aB/U^, and (2.50), we can cast (2.49) in the form 

r m C 

i -”y • (me + 85) 

P' ~ e G(q) (2.51) 


where C f is the skin-friction coefficient. 


1 c = T 

2 f TT 2 

P U 


Using the long-wave dispersion relation 


(2.52) 


w = Wq + 


and p' = p - p ro , we find that 
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(2*53) 


The quantity 1 1 /X q | is, of course, given to us by Figure 1. The skin 
friction coefficient could in principle be obtained from the actual 
computation at any time level. However, within our approximation one may give 
an estimate for Cf taken, for example, from boundary layer theory. Thus, 
for laminar flow and subsonic Mach numbers we may use the Blasius value of 

= 2/3/R^ . More on this in the section describing the numerical work 
which tested the efficacy of the boundary condition (2.53). 

We close this section by noting that the same methodology used here to 
derive outflow boundary conditions for the flat-plate problem can be used to 
extract boundary condition on the downstream side of other configurations such 
as the conditions in the wake of an airfoil, flow in a pipe, and so on. Also, 
other quantities have perturbations which satisfy (2.53), and one might try to 
apply the outflow condition to them as well. Another potential application is 
for unsteady flows where the dispersion relation for (2.13) - (2.15) would be 
carried out around some central wavenumber 6 C rather than around 8=0. 


3. NUMERICAL RESULTS 

In this section we present some numerical results for the boundary 
condition described above. These results are preliminary; more complete 
results will be presented elsewhere. Using the non-dimensionalization 
described above and the Blasius expression for Cf, the boundary condition 
(2.53) can be written as 
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fe " A(p - P») + B || (3.1) 

where 

A = " 9 1 x Q | ; B = “r 

We implemented (3,1) in a Navier-Stokes program using a version of the 
MacCormack scheme that is second-order accurate in time and fourth-order 
accurate in space [7], The numerical scheme was time consistent and did not 
use any acceleration technique such as a local time step. 

For our numerical experiment we used an = 0.4 boundary layer. The 
unit Reynolds number (per foot) was 3.0 x 10^. The inflow was taken at 1.0 
feet from the leading edge. We performed calculations using two different 
locations of the outflow boundary; x = 2.067 feet and 3.0 feet from the 
leading edge respectively. In both cases we used a grid of 31 points in the y 
direction (exponentially stretched). In the x direction we used 17 and 31 
grid points respectively for a mesh size of .067 feet. The inflow and initial 
data were obtained from a computer program which solved the boundary layer 
equations. Convergence to steady state was assumed when the maximum over the 
grid of \^\ was less than 10"^. 

The boundary condition (3.1) was compared with the condition (1.1). 
Equation (3.1) has the property that the steady state pressure is not forced 
to be equal to p. Thus, there may be some discrepancies in the steady state 
solution. The numerical experiments were designed to assess the behavior of 
(3.1) with regards to: 
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a) Convergence rate 

b) Robustness with respect to changes in the coefficients 

c) Quality of the steady state. 

We summarize our results below. 

Convergence to steady state was significantly accelerated by the use of 

(3.1) . The boundary condition (1.1) requires specification of a parameter a. 
The convergence rate is sensitive to the choice of a. In [5] an optimal a 
for a class of problems is presented. Because of the length of the computer 
runs, we were unable to determine the optimal a for the present problem. 
For the case of the outflow boundary at x = 3.0 feet from the leading edge a 
steady state was obtained using (3.1) in 30880 time steps. Using (1.1) and 

a = 1.45 the steady state required 97020 time steps. Using a = .6 which 
would be obtained from the formula in [6] based only on the free stream 
velocities and sound speed significantly degraded the convergence rate. For 
the case of the boundary at x = 2.067 feet from the leading edge, the use of 

(3.1) required 20960 time steps for convergence. The same run with a = 2.71 
required 52080 time steps. 

Robustness of the boundary condition (3.1) was tested by modifying the 
coefficient A in (3.1). A complete study has not yet been made, but 
decreasing A by a factor of 2 (outflow boundary at 2.067 feet) increased the 
number of time steps required for convergence by 560 and made only small 
changes in the final steady state. 

The quality of the steady state was assessed both by comparing solutions 
generated by the two boundary conditions and by comparing the solutions 
generated by the two positions of the outflow boundary. We first consider the 
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case of the outflow boundary at 3.0 feet. In this case the outflow pressure 
obtained from using (3.2) differed from by less than .04%. In Table I we 
indicate the value of U and V (horizontal and vertical velocities non- 
dimensionalized by U ) for the steady states obtained from (3.1) and (1.1). 
The data are presented at different x locations indexed by the grid point. 
(The outflow was at grid point 31.) In all cases the velocities are shown for 
the first grid point away from the wall as this was where the maximum relative 
difference occurred. The relative differences decreased significantly for 
grid points away from the wall where the velocities increased in magnitude. 

Using (1.1), we observed a slight oscillation in V near the outflow. 
This may be due to the fact that the outflow pressure is fixed at the free 
stream pressure while the flow appears to undergo a very small acceleration. 
For the boundary condition (3.1), V decreases rapidly close to the outflow. 
Thus, V near the outflow appears to be incorrect with both boundary 
conditions. Away from the outflow it is clear that V is the quantity that 
is most sensitive to the boundary condition although the differences are small 
and may also be affected by truncation error. The horizontal velocities are 
very close. At grid point 27 the two values of U differ by about 1.7%. 

In Table II we list the values of U obtained from using (1.1) and (3.1) 
with the outflow boundary at 2.067 feet from the leading edge (grid point 
17). These are compared with the values obtained when the outflow boundary 
was at 3.0 feet from the leading edge. We assume these values are accurate 
for the problem of the plate extending to infinity. Ideally, one would want 
the solution to be independent of the position of the outflow boundary. The 
values in these tables are the first grid point away from the wall where the 
relative differences are greatest. 



It is apparent from Table II that the values for U are very insensitive 
to the position of the outflow boundary* and to the choice of (1.1) or (3.1). 
The values obtained from (3.1) are slightly less accurate, but the differences 
are very small. At grid point 15 the relative error is only about 1.5%. 
These errors decrease for grid points away from the wall. 

The results in Table III show that there are significant boundary errors 
in V for both boundary conditions; however, these values are very small 
compared to the free stream values. At grid point 13 the relative error in V 
is about 8.5% for (3.1) and about 10.9% for (1.1). It would require a more 
complicated problem to assess the meaning of these differences. 

The boundary condition (3.1) has also been used to compute the steady 
flow field over a curved surface. A particular example is the flow over a 
curved hump with the curvature chosen so as to accelerate the flow from 

“ .7 to = 0.76. The height of the hump was .04 feet and the 
curvature extended over a range of 1.2 ft. The surface then became flat and 
the computational domain extended a further 1.2 ft downstream. The unit 
Reynolds number was 3.0 x 10 while the boundary layer thickness at inflow was 
.0095 ft. The Reynolds number based on the boundary layer height at inflow 
was 896.00. Using (1.1), and several different values of cx, convergence to 
the tolerance described above could not be obtained. Using (3.1) convergence 
in 52,760 iterations was obtained. 

Finally, we point out that the methodology described above can be applied 
to schemes that are inconsistent in time. As one example, we can consider 


systems of the form 
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where the flux functions F and (1 are the same as in (2.1) and E(U) is a 
preconditioning matrix (see [9] for appropriate examples). The matrix E(U) 
can be included in the linearization with only slight changes in the algebra. 

au„ - 

(Note that zero so that the only modification is that ^ is 

replaced by E(U Q ) |^-) . 

A simpler example is the use of a uniform grid in the viscous region and 
a different (larger) uniform grid in the inviscid region with the local time 
step appropriate for stability used in each region. It is easy to see that on 
the differential equation level this can be modelled by a matrix E which is 
I in the viscous region and of the form •— I where B is some large number 
depending on boundary layer thickness. This only affects the inviscid 
matching conditions ((2.22) - (2.24)). We have found that the parameters Xq 
and Uj are insensitive to B, if B = 0(e ^). 

4. CONCLUSION 

We have introduced a methodology which enables us to derive an outflow 
boundary condition for a range of steady flows. Preliminary results indicate 
that the resulting boundary condition may be effective in accelerating the 
computation of steady subsonic flows with little loss of accuracy and no free 
parameters to be determined. The method is also applicable to incompressible 
and supersonic flows. Extensions to unsteady flows are also possible. 
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Table I. Non-dimensionalized U and V obtained from (1.1) and (3.1) 
(y = .00088 ft.) 


X 

grid point 

U/U^ (1.1) 

U/U m (3.1) 

v/v ra (1.1) 

(V/u) w (3.1) 

ii 

.1234 

.1234 

.331 x 10“ 4 

.331 x 10~ 4 

15 

.1144 

.1144 

.263 x 10" 4 

.261 x 10" 4 

19 

.1071 

.1072 

.218 x 10“ 4 

.213 x 10" 4 

23 

.1009 

.1014 

.194 x io“ 4 

.174 x 10" 4 

27 

.09526 

.09680 

.201 x 10" 4 

.131 x 10" 4 

31 

.08978 

.09429 

.231 x 10" 4 

.0821 x 10" 4 


Table II. Non-dimensionalized U obtained from (1.1) and (3.1) with 

outflow boundary at x = 2.07 ft. compared with values obtained 
with outflow boundary 3 ft. from leading edge (y = .00088 ft). 


X 

grid point 

U/U^ (1.1) 

U/U (3.1) 

oo 

U/U^ outflow at 3 ft* 
from leading edge 

5 

.1420 

.1421 

.1421 

7 

.1348 

.1350 

.1349 

9 

.1286 

.1289 

.1287 

11 

.1232 

.1238 

.1234 

13 

.1184 

.1195 

.1186 

15 

.1139 

.1161 

.1144 

17 

.1106 

.1143 

.1106 
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Table III. Non-dimensionalized V obtained from (1.1) and (3.1) with 

outflow boundary at x - 2.07 ft. compared with values obtained 
with outflow boundary 3 ft. from leading edge (y = .00088 ft). 


X 

V/V* (1.1) 

v/v^ (3.1) 

V/V outflow at 3 ft. 

grid point 



from leading edge 




Figure 1 
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